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The Chemical Evolution of Globular Clusters - I. Reactive 
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ABSTRACT 

We propose a new chemical evolution model aimed at explaining the chemical prop- 
erties of globular clusters (GC) stars. Our model depends upon the existence of (i) 
a peculiar pre-enrichment phase in the GC's parent galaxy associated with very low- 
metallicity Type II supernovae (SNe II), and (ii) localized inhomogeneous enrichment 
from a single Type la supernova (SNe la) and intermediate-mass (4-7 M©) asymptotic 
giant branch (AGB) field stars. GC formation is then assumed to take place within this 
chemically-peculiar region. Thus, in our model the first low-mass GC stars to form are 
those with peculiar abundances (i.e., 0-depleted and Na-enhanced) while "normal" 
stars (i.e., 0-rich and Na-depleted) are formed in a second stage when self-pollution 
from SNe II occurs and the peculiar pollution from the previous phase is dispersed. 
In this study, we focus on three different GCs: NGCs 6752, NGC 6205 (M 13) and 
NGC 2808. We demonstrate that, within this framework, a model can be constructed 
which is consistent with (i) the elemental abundance anti-correlations, (ii) isotopic 
abundance patterns, and (iii) the extreme [0/Fe] values observed in NGC 2808 and 
M 13, without violating the global constraints of approximately unimodal [Fe/H] and 
C+N+0. 

Key words: nuclear reactions, nucleosynthesis, abundances - stars: abundances - 
stars: AGB and post-AGB - stars: chemically peculiar - globular clusters: individual: 
NGC 6752, NGC 6205, NGC 2808. 



1 INTRODUCTION 

As laboratories for both stellar and galactic physics, globu- 
lar clusters (GCs) rank among the valuable tools available to 
astronomers. Ancient, co-eval, equidistant, and apparently 
mono-metallic, GCs give the appearance of being elegantly 
simple aggregates of stars. Behind this "mask of simplicity" 
though, lurks perplexing observations which have defied ex- 
planation for several decades. 

First, each well-studied cluster to date shows star- 
to-star abundance variations of the light elements C, N, 
O, Na, Mg, and Al (JKraft 1994; Gratton ct al. 2004 and 
references therei n), not shared by field stars with simi- 
lar metallicities (jCratton et al.|[2000l ). These variations fol- 



over, ISmith et al.l (|2005D noted that stars in one of the 
Galactic globular clusters, NGC 6121, showed evidenc e for 
a Na-F anti-cor r elation , whereas iPasauini etall l|2005l ) and 
iBonifacio et~al] (|2007h found a Na-Li anti-correlation in 
NGC 6752 and 47 Tucanae, respectively. The abundances 
of Si, Ca, the iron-peak (e.g., Fe, Ni, and Cu) and neutron- 
capture elements (e.g., Ba, La, Eu) do not show the same 
star-to-star scatter, nor do they vary with the light elements 
(fCr atton eralll2004l : [james et al.ll2004l : lYong et al.|[2006bl . 
200d). The most important exception to this is the most 
massive cluster lj Centauri, that shows both, a spread in 
age and [Fe/H], and a rise in s-process element abundances 
with increasing [Fe/H]. Evidence fo r u) Cen suggests it 
may have a n extragalactic origin (e.g. Maiewski et al. 200' 



anti-correlated (e.g 'Kraft et al. 


'l993VGrundahI et al."200^; 
2003: Snedcn ct al. 200i: 


Cohen et all |2002; Yong et al 


YonE et all l2005l: Cohen & Melended 20051: ICarretta et al.l 


20061: iGratton et al.ll2007l:lMarino et al.ll2008l). ImDortantlv 



Smith et al. 2000: Gncdin et al. 2002; Bekki fc Frcema: 
20031 : iRomano et al.H2007l : iMarcofini et al.ll2007h . 



the sum of C-l-N-l-O is observed to be approximat ely con- 
stant (e.g. Ilvans et al] Il999l : ICarretta et al] l2005f ). More- 



The above abundance trends have been found for 
stars in all evolutionary phases, from the main-sequence 
turn-off through to the tip of the first giant branch, 
which has lent support to the "self-pollution" hypothesis 
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l|Cottrell fc Da Costa|[l98ll : lDantona et al.lll983l ). This is in 
contrast to deep mixing, where the abundance anomalies are 
produced by internal mixing during the ascent of the giant 
branch, after the f i rst dredge-up (,Sweigart fc Mengel. 19791 : 
|Pinsonneaui3ll997l : ICharbonneJll994 ). Further evidence for 
self-pollution comes from a lack of variation in the O-Na and 
Mg-Al anti-correlations with luminosity as the stars ascend 
the first giant branch. 

According to the self-pollution scenario a previous 
generation of stars "contaminated" the atmospheres of 
stars observed today in GCs, or provided much of the 
mate r ial from which those s t ars formed (e.g. Jehin et al.l 



1 19981 : iParmentier et al.l 1 19991 : iTsuiimoto et al.1 



20071 ). The 



roughly constant [Fe/H] abundances in a given GC led 
to the assumption that the polluters were low-metallicity, 
intermediate-mass asymptotic giant br anch (AGB) stars, 
with masses between ~ 4 and 7 M© (|D'Antona fc Caloil 
12004 iBekki et al] l2007h . The hot bottom burning ex- 
perienced by these stars provides the hydrogen burning 
environment (at least qualitatively) that can alter the 
abundances of the light elements. Detailed AGB mod- 
els have most ly failed to explain the observed abun- 
danc e trends (iDenissenkov fc Herwid [20031: iFenner et all 



l2004l : ICohen fc Melended l2005l : iKarakas et al.1 l2006al ). 
but model uncertainties, including convection and mass 
loss, render predictions uncertain (jVentura fc D'Antonal 
l2005al 3). Recently, slow- winds from rotating massive stars 
(IPrantzos fc Charbonnelll2006l : ISmithll2006l : iDecressin et~all 
l2007l ) have also been proposed as another possible source of 
the abundance anomalies. 

The canonical AGB enrichment scenario consists of two 
stages in which the first stars to form are those pre-poUuted 
by core collapse supernovae (SNe II), and have [0/Fe] ~ 0.5 
and [Na/Fe] ~ —0.1. After the first SNe II have driven away 
the gas from the proto-cluster, gas polluted by AGB winds 
starts to accumulate at the centre of the cluster and a sec- 
ond generation of stars with peculiar abundances can form. 
These second- generation stars will have depleted F and O, 
along with enhanced in N and Na with respect to the first 
generation. Variations to this scenario include (1) the AGB 
ejecta-enriched gas is al so accreted on to the su rface of some 
(e.g., segregated) stars jParmentier et al]|l999l ). and (2) the 
AGB pollution comes, not only from the first generation of 
stars, but also from all AGB s tars of the sa tellite galaxy 
where the GC is forming (Bekki et al.1 bOOTl ). The second 
assumption was made to allow for a larger reservoir of ma- 
terial for the formati on of the second generation of stars 
l|Bekki fc Norrisll2006h . as if all AGB pollution comes from 
the first generation of stars, is very difficult to explain the 
high fractio n (as high as 50%) of peculiar-to- normal stars 
in GCs (e.g..lD'Antona fc Caloill2004l:lD'Antona et al.ll2005l : 



ICarretta et al.ll2006l : iPiotto et al.ll2007l ). The same results 
can be achieved by assuming that the GCs were much more 
massive in the past (by a factor of ~10-100x) and lost pref- 
erentially most (90-99%) of their first generation of stars 
through tidal intera ctions with the Milky Way (MW; e.g. 
iD'Ercole etHlbOOSl ). 

Another way to solve the problems associated with 
the AGB-"pre-pollution" scenario is to assume that the 
polluters were not A GB stars but Type la supernovae. 
iMarcolini et al.l (|2006l ) suggested that the O-depleted stars 
observed in dwarf spheroidal galaxies (dSphs) as well as in 



the peculiar globular cluster uj Gen l|Marcolini et al.|[2007f ) 
can be explained if those stars were born in a small iron-rich 
volume inhomogeneously polluted by a single SN la. The 
ejecta of SNe la is very rich in Fe and extremely poor in O 
(compared to SNe II) thus stars forming in such confined 
regions would show very low values of [O/Fe] (usually up 
to ^ —0.5). Dwarf galaxies, however, show a large spread 
(1.5 dex) of [Fe/H], hence it is not clear if this model is 
applicable to the problem of the abundance anomalies in 
mono-metallic GCs. 

In this paper we present a GG chemical evolution 
model where all of the observed chemical peculiarities arise 
from the inhomogeneous pollution of a single SN la and 
intermediate-mass AGB stars. This pollution is mixed with 
gas that was pre-enriched by very low-metallicity SNe II. 
Since each of these so-called "polluters" precede the forma- 
tion of the GC itself, we call this type of enrichment "ex- 
ternal pollution" (according to Bekki et al. 2007) from the 
Milky Way (or host galaxy). When the first stars with pe- 
culiar abundances are forming, SNe II start to explode, self- 
polluting the ISM and dispersing the inhomogeneous region. 
As star formation proceeds, the abundances of the forming 
stars will move from "peculiar" to "normal", where normal 
in this context means stars with chemical properties typical 
of SNe II enrichment. Eventually, the SNe II explosions will 
cease any further star formation (SF) and the GG will evolve 
passively. In previous studies, SN la had been discarded as 
polluters because the [Fe/H] is essentially constant within a 
given GC, and SN la are efficient producers of iron. Never- 
theless, we will show that our model is able to reproduce the 
observed anti-correlations, while keeping the iron content of 
the newly forming stars approximately constant. 

In the following we test our scenario on a GG with a typ- 
ical [0/Fe]-depletion pattern (on the order of ~ —0.5), simi- 
lar to to the well-studied cluster NGG 6752 ([Fe/H]=-1.56). 
In addition we test our scenario on a more extreme case 
such as that of NGC 2808 ([Fe/H] =-1.15) and M 13 
([Fe/H] =-1.50) jHarrid [199^ ). with [O/Fe] -depletions as 
large as 1 dex. 



2 INITIAL CONDITIONS 

It is generally assumed that GCs were assembled during the 
formation of the Milky Way halo (or their parent galaxy) 
even though their chemical properties differ from those halo 
field stars of the sam e metallicity (e.g.. ICratton et al.ll2000l : 
iFreeman fc Bland-Ha wthorn 200^). It is likely that part, or 
all, of this difference can be traced to the very different con- 
ditions experienced by the respective systems e.g., deeper 
potential wells and the associated burst of SF. 

Following this framework, we assume that a typical GC 
forms initially within gas pre-enriched by field polluters such 
as SNe II, but also from SNe la and AGB stars associated 
with (parent) halo star formation. However, due to the early 
stage of the formation of the halo (with respect to the rest of 
the parent /Milky Way), SNe la and AGB stars are not im- 
portant in shaping the mean chemical properties of the gas 
unless they are polluting it inhomogeneously (i.e., polluting 
a small localized region) . Thus, we assume that a typical GC 
will form initially in a peculiar region of the halo that was 
pre-poUuted by SNe II, and in which the inhomogeneous pol- 
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lution of a single SN la and intermediate-mass AGB stars is 
added. In our model, all the chemical peculiarities observed 
in GC stars (i.e., O-depleted and Na-rich stars) which are 
not shared by their field counterparts arise due to this pe- 
culiar localized inhomogeneous pollution. 

A schematic representation of how these initial con- 
ditions can be achieved is shown in Fig [T] According to 
our model, AGB material is deposited in the vicinity of a 
SNe la progenitor binary system. As the secondary explodes 
as a SN la and the supernovae remnant (SNR) expands, it 
will collect this surrounding AGB material in its shell. The 
numerical e volution of a radia tive SNR has been analyze d 
in detail bv lCioffl et al] (|l988l ) and iThornton et al.l (|l998D . 
These studies found that for normal ambient gas density of 
no = 1 — 10 cm~"^, the evolution of a SNR lasts for a few 
Myrs and that it reaches radius of ~50-90 pc (with the larger 
radii achieved for the smaller ambient gas densities). After 
radiating away most of its energy, the remnant will first stall, 
lose its identity and then collapse back (e.g. ISlavin fc Coxl 
1 19921 . ; approximately after Rmax/cs ~ 5-9 Myr, assuming 
an ambient sound speed Cs of 10 km s~^). The result of 
this evolution will be a localized region of the halo, polluted 
inhomogeneously with SN la and AGB ejecta. These will 
be the initial conditions of our model which we define as 
the "pre-peculiar enrichment" scenario; we assume that the 
proto-GC starts forming stars in this peculiar region. 

Regions polluted inhomogeneously by SNe la arise nat- 
urally in the hydrodynamica l simulations of dSphs chemi- 
cal enrichment performed by iMarcolini et al.l (|2006l . |2008| ) . 
ICescuttil (|2008l ) also showed that a model which takes into 
account inhomogeneous pollution is able to explain the ob- 
served spread in s- and r-process ele ments at low metal- 
licities in the Galactic halo (see also llshimaru fc Wanaiol 
ll999l:lArgast et al!'200(f). Since both the peak SN e la rate 
I Matteucci fc Rccchi 2001; Ma nnucci et al.l \200dt) and the 
lifetime of a low-metallicity 5 M p, star (e.g., Schaller et al.l 
1 19921 : iKarakas fc Lattanzlol l2007l l are comparable (~ 80 — 
100 Myr), this apparently pecu l iar e nrichment scenario 
seems reasonable. iMannucci et al.l l|2006l ) proposed that half 
of the SNe la associated with a starburst should explode 
on a timescale of the order ~ 100 Myr. Marcolini et al. 
(2006,2008) have shown that even if a timescale greater than 
1-2 Gyr is needed to appreciate the cumulative effect of 
SNe la in the mean [a/Fe] ratio of spheroidal galaxies, the 
effect of a single SN la inhomogeneously polluting the ISM 
could be important i n the first few hundred Myrs (see also 
IMarcolini et al.ll2007l . for and application of this scenario to 
the particular case of w Gen). 



3 TWO-REGION MODEL 

Our model is based on what could be termed a "two-region" 
chemical evolution framework with a simplified prescription 
for the chemical enrichment. As a first attempt to reproduce 
the main chemical properties of GC stars, we do not follow 
the ejecta of different mass progenitors, but we use IMF- 
averaged yields for both intermediate-mass AGB stars and 
SNe ejecta (see |4] for more detail). A sketch of the model 
is shown in Fig. O where the two regions - "inner" and 
"external" - are shown. The dark blue "inner region" in 
Fig. |2] corresponds initially to the blue-green circle in the 



^ c 
^ o o o 






Figure 1. Schematic framework leading to the initial conditions 
required for our model: 1) low-metallicity 4—7 Mq AGB stars 
(green dots) start polluting a particular region of the halo which 
has been pre-enriched by SNe II; 2) a SN la explodes near this 
region and during its expansion collects the AGB material in its 
shell; 3) the SN remnant radiates away its energy and the gas 
starts re-collapsing, mixing the SN la ejecta with material from 
the AGB stars; 4) the final result is a peculiar region in the Galac- 
tic halo, where on top of the mean pre-enrichment from early 
SNe II, the inhomogeneous contribution by a single SN la and a 
number of intermediate-mass AGB stars is added. These are our 
assumed initial conditions for proto-GC formation. Note that the 
light-blue background volume (mainly polluted by early-SNe II) 
represents the mean chemical properties of the gas out of which 
"normal" field halo stars are formed. 



last panel of Fig. [T] polluted inhomogeneously by a SN la 
and intermediate-mass AGB material; we will assume that 
the SF related to the GC formation only takes place inside 
this region. In the following, we will use the term "inho- 
mogeneous pollution" to denote that this region has differ- 
ent chemical properties from the surrounding medium (note 
that the chemical properties inside the inner region itself are 
homogeneous) . 

The external region is simply described by the initial 
value of [Fe/H]isM given in Table [1] which should refiect 
the mean chemical properties of the proto-halo gas at the 
epoch of GC formation (enriched mainly by SNe II). The 
inner region is initially described by its radius (Rin i.e., the 
degree to which the SN la was initially confined) and the 
mean chemical properties within, which is set by the num- 
ber of 4—7 Mq AGB stars, A^agb. Since most of the iron in 
the inner region is provided by a single SN la the value of 
Rin also defines the value of [Fe/H]in in the inner region. It 
is obvious that the localized iron-rich ejecta from the SN la 
will substantially increase the value of [Fe/H] inside the in- 
ner region compared to the value in the external region, 
whereas intermediate-mass AGB stars increase the N, Na, 
the neutron-rich Mg isotopes and, possibly, Al. A glance at 
Table [1] shows that the inner region is a factor of ~ 10— 100 
more [Fe/H] rich than the external ISM and also has a sim- 
ilar higher overall metallicity. For this reason the first gen- 
eration of stars will born with peculiar chemical properties: 
they will be depleted in O and Mg (because of the single 
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Table 1. Initial conditions of the model. Note that only three 
of the parameters are free as the value of [Fe/H]in is defined 
by (see text for more details). 



Model 


[Fc/H]isM 


[Fc/H]i„ 


Rin (pc) 




NGC 6752 


-2.55 


-1.60 


36 


250 


M 13 


-3.50 


-1.50 


31 


170 


NGC 2808 


-2.85 


-1.10 


24 


180 



SN la) and they will be enhanced in N, Na, and Al (due to 
AGB polluti on). Note that in c anonical AGB self-pollution 
models (e.g. iFenner et al.l |2004| ) these stars are the last to 
form. While a detailed discussion about each element will 
be presented in the following sections, here we would like to 
point out that, in general, the trend for SN la is to decrease 
the [Xi/Fe] value when Xi is not (or negligibly) produced in 
intermediate-mass AGB stars (e.g., O and the elemental Mg 
abundance to a lesser extent). In the case when the elements 
are overproduced in AGB stars (e.g., N, Na and, possibly, 
Al), the effect of the SN la is to mitigate their enhancement. 

We assume that the star formation takes place only in- 
side the inner region. As the SF proceeds, massive stars 
(M© ^ 8 Mq) start exploding (as early as ^ 3 — 4 M yr 
after the onset of star formation ^Scha ller et al.l Il992l ) as 
SNe II, ejecting freshly synthesized metals. We assume that 
after each SN II explosion the inner region expands, driven 
by the explosions themselves, mixing the inner-region gas 
with the surrounding ISM, while the SF is still proceedings 
inside this region. During the evolution, we simply assume 
that the density remains constant and thus the mass inside 
a fixed radius is M= |7rpoR^. For each of the proposed GC 
models we assume a density of po = 4 x 10~^* g cm""^ which, 
for example, corresponds to a total mass of Mism — 3 x 10^ 
M© inside a 500 pc radius sphere. Even if we were to as- 
sume a lower SF efficiency, there would still be more than 
sufficient gas to form a typical GC. 

The simple relations governing the evolution after each 
SN II explosion can be summarized as: 

R{t + At) = R{t) X /exp (1) 

X,{t + At)^X,{t)+X,{AR) + {X,{SNII)) (2) 

where /oxp is a constant governing the logarithmic expansion 
of the inner region after the explosion of every single SN II, 
At is the time interval between two successive SNe II explo- 
sions, Xi{t) is the mass fraction of an element X at a time 
t within the radius R(t), Xi(AR) is the mass fraction of an 
element X initially present in the corona R(t-)- 1)— R(t), and 
finally, {Xi{SNII)) is the average mass fraction of element 
Xi ejected by a single SN II (see Sec. |4]). For example, the 
evolution of the [Fe/H] abundance in the inner region at a 
time {t + At) (where At is the time between two successive 
SNe II) is computed by adding the Fe initially present in 
this region at R(t) to the iron ejected by the newly explod- 
ing SNe II and to the iron initially present in the external 
region R(t -I- At)— R(t). This number is then divided by the 
total amount of hydrogen inside the new radius R(t -I- At). 
All the other variables are calculated accordingly. We as- 
sume that this expansion is described by a logarithm law 
with a free parameter (usually of the order of /oxp ~ 1.001) 



normal pollution 
(mainly SNe II) 




Figure 2. Schematic of the model. Initially, a localised volume 
(inner blue region) arises with the pollution from a single SN la 
and multiple intermediate-mass AGB stars. After a new genera- 
tion of stars is born the associated SNe II begin to pollute and 
expand the inner volume, as well as mix it with the surrounding 
lower- [Fe/H] interstellar medium. 



which fixes the increase of R(t) after each SN II explosion 
(i.e., R(t + At) =R(t) X /exp). 

In Fig. [2] the dashed lines represent a sketch of the evo- 
lution of the inner region and the fact that the initial con- 
tribution of the single SN la (as well as AGB stars) is di- 
luted during its expansion and the subsequent mixing with 
the outer, lower metallicity ISM. As soon as the first SNe 
II start to explode there are two contrasting mechanisms 
shaping the chemical properties of the newly forming stars. 
First, the freshly ejected SNe II metals directly pollute the 
gas where new stars are forming. Second, the inner volume 
is expanding and the gas inside this region is mixed with 
lower metallicity material. In Sec. 15.11 we will show that, 
during the chemical evolution process, it is possible to keep 
the [Fe/H] abundance constant in the inner region (where 
the new stars are forming) owing to these two contrasting 
mechanisms. 

While stars with peculiar properties form during the 
first stage, during the second stage the chemical properties 
of the forming stars a -enhanced) as SNe II start exploding. 
iMarcolini et al.l l|2006l ) showed, for the specific case of dwarf 
galaxies, that even if a single SN la can temporarily pollute 
a small volume and dramatically change its chemical com- 
position, its imprint disappears almost completely once the 
gas inside this volume mixes with the ISM. In the next sec- 
tion we will show that the last stars forming in our model 
(stars with [0/Fe]~ 0.5) do not show any signature of the 
former AGB or SN la pollution, having chemical properties 
typical of pure SNe II ejecta. 

We will also show, in Sec. 15.11 that in our model the 
inner region is more metal rich than the outer region and, 
accordingly, the cooling is much more efficient there. There- 
fore, this region is more likely to host an extremely enhanced 
SF - typical of GCs - than the outer ISM (where normal 
halo stars are forming). 

Numerous authors have discussed, from an energet- 
ics perspective, the possibility that stellar clusters can 
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Table 2. Mean SNe II stellar yields averag ed over the progenitor mass range 10—60 Mq for a Salpeter I MF a nd for different 
authors: Wfcw jWoosIev fc WeaveJ l|l995i ): CfcL= IChieffi fc Lim ongi (2004); KOB= lKobavashi et al.l l|2006l 'l. The yields of 
Fe and He are given in solar masses while for different elements we show the [Xj/Fe] ratios. 



SNe II model 


Fc 


He 


[C/Fc] 


[N/Fc] 


[O/Fc] 


[Na/Fc] 


[Mg/Fe] 


[Al/Fe] 


W&W (Z=0.0002 ) 


6.10-2 


6.7 


0.13 


-1.73 


0.51 


-0.39 


0.15 


-0.24 


W&W (Z=0.002) 


6.9e-2 


6.6 


0.07 


-0.86 


0.46 


-0.16 


0.16 


-0.04 


C&L (Z=0.0001) 


l.Oe-1 


5.7 


0.41 


-1.91 


0.63 


0.05 


0.61 


-0.14 


C&L (Z=0.001) 


l.Oe-1 


5.7 


0.40 


-0.96 


0.63 


0.26 


0.61 


0.11 


KOB {Z=0.001) 


7.5e-2 


5.7 


-0.16 


-0.80 


0.64 


0.05 


0.63 


0.15 


KOB {Z=0.001+HN) 


9.3e-2 


5.7 


-0.26 


-0.89 


0.52 


-0.08 


0.52 


0.28 


Model (SNo II) 


9.0O-2 


6.0 


-0.10 


-0.90 


0.50 


-0.20 


0.50 


-0.10 



experience SNe II self-enrichment, and there is a wide 
ranee of evidence as to whether, and to what ex- 



tent, this is possible in GCs (e.g. 


Dopita & Smith' 'l986l; 


Morgan & Lake 1989; SmithI 19961; 


Parmontior & Gilmore 


2001 


; iGnedin et all 12002 




Thoul et all I2OO2I; IParmentier 


2004 


; Recchi & Danziger 


20051; IPrantzos & Charbonnel 


2ooe 


). For example, it has been suggested that stellar winds 



and SNe explosions from the first generation of stars can 
form a metal-enr iched supershell, where f urther star forma- 



tion is triggered llBrown et al.|[l99ll.ll995l;lThoul et al.ll2002l; 
iParmentieij |2004| ; iRecchi fc Danziged |2005| ). iBrown et all 



119911) 

concluded that a small entirely self-enriched sys- 
tem of stars typically achieves a metallicity in the range 

10~^ 10~^ Zq, and that the second generation of stars 

would be expected to be extremely homogene ous in com- 
position. iMelioli fc de Gouveia Dal Pinol l|2004r ) studied in 
detail the heating efficiency (i.e., the fraction of SNe energy 
which is not radiated away) of SNe exploding in a starburst 
environment with properties similar to that of our model. 
These authors found that the heating efficiency can be very 
close to zero (i.e., all the SNe energy is efficiently radiated 
away) for, up to, a few tens of Myr, before rapidly increas- 
ing to unity and leading to a possible galactic wind. While 
we cannot follow the dynamical evolution (and the energet- 
ics) of our model, it is possible to speculate that this second 
stage of star formation could be as rapid as ~ 20 — 40 Myr 
(with SNe II directly polluting the newly forming stars), 
or it could proceed in a long series of weak bursts if the 
star formation is self-regulated - i.e., the gas is temporar- 
ily expelled out of the star forming regions by the SNe and 
re-collapses (due to the outer pressure) once the SNe II en- 
ergy is radiated away. If the star formation is self-regulated, 
the SF history can be as long as ~100 Myr. A longer SF 
history would require the contribution of AGE and SNe la 
self-pollution, which does not seem to be required in our 
model, to explain the chemical properties of GCs. Finally, 
in our case, the formation of a galactic wind would be associ- 
ated with the cessation of the SF itself and it would happen 
after ~ 1 — 5 x 10'' SNe II explosions, depending upon the 
details of the model. Although these earlier studies seems 
to support the formation of second generation of stars after 
the occurence of an explosive event, and the formation of 
a bound system, hydrodynamical simulations of the SNe II 
self-enrichment phase in quantitative detail are needed to 
further quantify the star-formation timescales and the via- 
bility of our model from a dynamical point of view. These 
will be presented in a forthcoming paper in this series. 



4 YIELDS 

As noted earlier we used IMF-averaged yields for both 
SNe II and intermediate AGB polluters. For SNe II, we 
explore the effect of using yie l ds computed by differ- 
ent a uthors |Wooslev fc Weaverl 1 19951 ; IChieffi fc LimongH 
12004 iKobavashi et all |2006| . ~e Table HJ and use the 
limits 10-60 Mp). Previous studies have found that the 
IWooslev fc Weaver! l| 19951 ) models tend to overproduce the 
amount of Fe - with a subsequent underestimation of 
the [a/Fe] abundances - for this re ason we half Fe val- 
ues as done in several studies (e.g., Timmes et al.l Il995l ; 
iGoswami fc PrantzosI I2OO0I ; iFenner et al.ll2004l ). While the 
set of yields are in reasonable agreement with one another 
for most of the elements , occasional differe nces as large as 
0.4 dex are present (e.g.. iGibson et al]|l997l ). We computed 
a set of yields that we use throughout the paper unless 
stated otherwise. These yields are the average value of the 
above mentioned theoretical set of yields when the agree- 
ment between them is good (when they difi'er less than 
0.4 dex). Otherwise, we chose the yields that best reproduce 
the low metalhcity (—2.0 ^[Fe/H]^ —1.0) halo stars (e.g 
iGratton et~alll2000l ). These computed values are labelled as 
"Model" yields in Table H 

For SNe la we use the yields of llwamoto et af] l|l999l ). 
specifically their case WDDl (note however that for the el- 
ements considered in this paper a different choice of their 
SNe la model would not affect the results). For a consis- 
tent chemical evolution model, the metallicity of the pre- 
poUuting SN la and AGB stars should be be roughly the 
same as the external region ([Fe/H]isM), while the metallic- 
ity of the self-polluting SNe II match that of the inner region 
([Fe/H ]in). Unfortunately, the set of yields of llwamoto et all 
( |l999l ) are only available for solar metallicity and we are 
not aware of low-metallicity SNe la yields available in the 
literature. 

In Table |3] we summarize the mean yields for 
intermediate-mass AGB s tars, calculated by av e raging the 
Z = 0.0001 values from iKarakas fc Lattanzij (|20071 ) us- 
ing a Salpeter (1955) IMF in the mass range 4-7 Mq (we 
will refer to this model as the "reference model"). As a 
comparison we also show the yield values obtained aver- 
agi ng similar datasets at the sam e metallicity published 



ag mg simuar datasets at the sam e metallicity publisned 
bv iKarakas fc Lattanziol (120031 ) an d llzzard et all 120071) . as 
well as the yields calculated by iFenner et al.l (|2004l ) for 
the specific case of AGB self-pollution in NGC 6752 at 
[Fe/H]= —1.40. In general, each set is in good agreement 
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Table 3. Mean SNe II, SNe la and intermedi ate-mass AGB y ields in solar masses. SNe II yields refer to the model in Table[2] 
SNe la yields refer to the WDDl model of llwamoto et alj |l999i). The r nean AGB yields are averag ed over the progeni- 
tor m ass range 4-7 Mp) u sing a Salpctor IMF for different authors: K A07= lKarakas fc Lattanziol ||2007|) : FE04= lFenner et alJ 
||2004|): KA03=jKarakas fc Lattanzio (200i); IZ07= llzzard et alJ l|2003); all the yields are for Z = 0.0001 except for the case of 
iFenner et al. which is at [Fc/H]=— 1.40. The yields of Fe and He are given in solar masses. 



Polluter 


Fe 


He 


C 


N 


O 


Na 


Mg 


Al 


Model (SNe II) 


9.00C-2 


6.00 


1.59C-1 


7.66C-3 


1.82 


1.60C-3 


1.49e-l 


3.23C-3 


Model (SNe la) 


6.72e-l 


5.660-3" 


5.42e-3 


2.85e-4 


8.82C-2 


8.77e-5 


7.69e-3 


4.38e-4 


AGB (KA07) 


2.27e-5 


1.45 


6.21e-3 


6.89e-2 


1.44e-3 


1.60e-3 


2.17e-3 


7.01e-5 


AGB (FE04) 


X 


X 


8.54e-3 


5.76e-2 


1.86e-3 


5.95e-4 


1.52e-3 


6.13e-5 


AGB (KA03) 


X 


X 


X 


X 


X 


7.79e-5 


1.70e-3 


1.04e-4 


AGB (IZ07) 


X 


X 


X 


X 


X 


2.92e-5 


2.15e-3 


2.02e-4 


Model (AGB) 


2.27e-5 


1.45 


1.55e-3 


6.89e-2 


1.44e-3 


4.00e-4 


2.17e-3 


3.50e-3 



" the value of He for SNe la is from iKobavashi et al ] l|2006h . 



although there are significant differences present for some 
elements (e.g., Na yields span ~ 2 orders of magnitude). 
A detailed study of how the production of the Ne, Na, 
Mg, and Al isotopes in AGB stars are effected by r eactio n 
rate uncertainties was carried out by llzzard et al.l |20o3). 
They concluded that the most uncertain yields are those of 
■^^Al and ^^Na, with variations of two orders of magnitude. 
The yields of '^''Mg and ^'^Al have typical uncertainties of 
one order of magnitude. These uncertainties are compara- 
ble to those resulting from stellar modelling uncertainties, 
where differences of one order of magnitude or more can 
be found by varying the mas s-loss r ate or convective mode l 
IVentura fc D'Antonal (|2005al lbl) and lKarakas et al.l (|2006bh . 

For this reason we also decided to use a set of ad — hoc 
intermediate-mass AGB yields that represent the different 
element production that we would require in order to re- 
produce the observational constraints in the framework of 
our model. This set of yields is presented in Table [3] as 
"Model". These values are within the limits of the above- 
mentioned set of yields, and always within a factor of four 
of the values of the reference model (with the exception 
of Al see Section 15.31). In the following, the models using 
iKarakas fc Lattanziol (|2007^ (reference model) will be shown 
as solid lines, while the models using our proposed set of 
yields will be shown as dashed lines. 

We would like to stress that the AGB yields we propose 
in Table [3] should not be taken as absolute values but as rel- 
ative ratios between the different elements. This is because 
it is possible to re-scale the yields by changing the num- 
ber of Nagb stars needed to fit the observations. We use 
nitrogen as a reference element and re-scale the other ele- 
ments accordingly. We choose nitrogen because it is better 
constrained than the other elements (e.g., Na, Al) since its 
production does not depend as strongly on reaction rate un- 
certainties. N ote that the recent revi sion of the ^'*N(p,7)^^0 
reaction rate (|Bemmerer et al ] |2006l ). ~ 60% slower at tem- 
peratures below lOOxlO^K, would result in even more N 
production. Nitrogen yields are, however, effected by other 
modelling uncertainties such as convec tion, mass loss and 
the modelling of th e third dredge-up (|Marigo et al.l l2003l : 
iKarakas et al.|[2008l ). For example, an increase in the mass- 
loss rate would result in smaller N yields because of a short- 
ened HBB phase (and also less Na, Mg and Al). 
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Figure 4. The predicted trend of [Fe/H] versus [O/Fe] for three 
models wi th different /exp v alues compared with the observational 
dataset of lYong et all ([2005) for the case of NGC 6752. The initial 
values of /exp are 1.0003 (dashed line), 1.0009 (solid line) and 
1.0027 (dotted line) and are increased to 1.0012, 1.0035, and 1.015 
after 750 SNe II have been exploded. The models are stopped after 
~1600 SNe II have been exploded polluting the ISM. It is evident 
that changing the value of /exp within a factor of three effects the 
[Fe/H] ratio of the forming stars by only 0.1-0.2 dex. 



5 RESULTS: THE CASE OF NGC 6752, M 13 
AND NGC 2808 

To test our framework in detail, we will focus on three 
weU studied GCs: NGC 6752, M 13, and NGC 2808. All 
these GCs have an intermediate metallicity ranging from 
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Figure 3. The predi cted trend of [Fe/H] versus [O/Fe] compared witii tiie observational dataset of lYong et al. I bOO^ (left panel) and 
ICarretta et al.l l[2003) (right panel) for the case of NGC 6752. It is possible to maintain an invariant [Fe/H] abundance owing to a 
tuning between the two opposing mechanisms: SNe H eject freshly synthesized Fe (and O), which is mixed with gas from the inner 
(inhomogeneous) volume along with material from the external ISM, which has a lower [Fe/H] value. 
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Fi gure 5. The predicted trend of of [Fe/H] versus [O/Fe], for the case of NGC 2808 (left panel), compared with the obse r vation al dataset 
of ICarretta et al.l (l2006[i . The case of M 13 (right panel) is plotted against the observational dataset of ISneden et al.l (|2004| , asterisks) 
and ICoheiT'fc'^eMndez 1I2OO5I. diamonds). 



-1.60 ^[Fe/H]=S -1.10 (e.g.,|Harry[l99^. While NGC 6752 
can be defined as a "normal" GC, NGC 2808 and M 13 show 
extreme values of the [O/Fe] ratio (down to ~-1.0) and are 
excellent cases to test our new proposed framework. In this 
section we show how different choices of the initial condi- 
tions in the pre-enrichment during the MW halo formation 
can led to different chemical pecuharities in the forming GC. 

5.1 [Fe/H] and Metal content 

Before going any further, it is important to follow the 
iron content of the forming stars in our model. As previ- 



ously mentioned, the possibihty of SNe II and SN la self- 
enrichment in GC have not been considered in the literature 
before, because they produce large amounts of Fe, whereas 
the [Fe/H] content in most of the GCs appears to be con- 
stant. 

The initial conditions of NGC 6752 are summarized in 
Table[T] the Rin of the inhomogeneous region polluted by the 
SN la is 36 pc and this volume is also poUuted by 250 inter- 
mediate 4-7 Mq AGB stars. The mean iron content of the 
external region is [Fe/H]= —2.55, reflecting the early stage 
of the halo formation, while the effect of the SN la increases 
the value of [Fe/H] up to —1.60 inside the inner region. In 
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Figure 6. Evolution of tiie metallicity \og{Z / Zq) versus [O/Fe] in the inner region for the three reference models: NGC 6752 (left 
panel), M 13 (central panel), and NGC 2808 (right panel). The dotted line corresponds to the case where we ignored the contribution 
of the SN la+AGB pollution only to the total metallicity log(Z/ZQ) (note that low initial [O/Fe] value can be achieved only due to the 
localized contribution of the SN la), while dashed lines correspond to the model with no pollution from AGB stars. Note that the initial 
value of \o^[Z / Zq) (thus at low values of [O/Fe]) for the dotted line corresponds to the value of the metallicity for the external region. 
Most of the metals in the inner region are initially produced by AGB stars while at the end of the evolution only SNe II contribute to 
the metal content. 



the same table, we also summarize the initial parameters for 
the cases of M 13 and NGC 2808. 

For each GC, the initial conditions are constrained by 
the observed stellar values of [Fe/H] and the peculiar values 
of [0/Fe]niin and [Na/Fe]max (or alternatively the value of 
[N/Fe]niax), as each of these values depends mainly on one 
polluter. Since most of the iron in the inner region is pro- 
duced by SN la, the initial [Fe/H] in content of the forming 
stars is set by the value of Rin, while oxygen is produced 
mainly by SNe II and Na by intermediate-mass AGB stars 
(see Tables [2] and |3} ■ We assume that a SN la always deposits 
the same amount of Fe (0.74 M©), hence decreasing (increas- 
ing) the value of Rin has the effect of decreasing (increasing) 
the hydrogen content inside the inner region and increasing 
(decreasesing) the [Fe/H] content (this is evident in Table [U 
where the value of [Fe/H]in increases for lower values of Rin). 
Once the [Fe/H] content in the inner region is set, the over- 
all [Fe/H]isM value constrains the initial [0/Fe]min of the 
model: this is because, according to a pure SNe II pollution 
model, the O content is proportional to the Fe content (and 
the Fe produced by SNe II can be ignored in the inner re- 
gion). At this point, the observational value of [Na/Fe]max 
is sufficient to constrain the number of intermediate-mass 
AGB stars needed to reproduce the Na-O anti-correlations. 

For example NGC 6752 and M 13 have similar initial 
values of [Fe/H], at ~ —1.60 and [Fe/H]~ —1.50, respec- 
tively, that reflects comparable values of Rin- In contrast 
the "external" values of [Fe/H]isM for each cluster differ 
by an order of magnitude (see Table [T] reflecting, as shown 
later, their different extreme values of [O/Fe]. In contrast, 
NGC 2808 shows an even higher value of [Fe/H]~ -1.10 be- 
cause of its smaller value for Rin=24 pc, even if its [Fe/H]isM 
is slightly lower than for the case of NGC 6752 (see Table [l|. 

In Fig. |3]we show the [Fe/H] versus [O/Fe] abundance 
evolution o f the r efere nce model for NGC 6 752. Data from 
lYong et al.) (|2005l ) and lCarretta et all (|2007l ) is included for 
comparison. To compare our model with the Carretta's data, 
we apply an offset of 0. 1 dex to our predicted value of [Fe/H] . 
This is to take into account a similar difference in the mean 
[Fe/H] values between the two datasets. It is clear from the 
figure that the [Fe/H] content remains relatively constant 



during the whole evolution and it is consistent with the ob- 
servational data. Likewise, the evolution of [O/Fe] increases 
monotonically from its initial peculiar value [0/Fe]~ —0.4 
to 0.5. 

The [Fe/H] abundance stays constant during the evo- 
lution because of the tuning between the two mechanisms 
changing the Fe content of the forming stars. While the 
newly exploding SNe II eject freshly synthesized Fe, increas- 
ing the [Fe/H] content, their action also has the effect of 
expanding the inhomogeneous inner region, mixing the in- 
ner material with the external ISM (with a lower [Fe/H] 
content). Both these two contrasting effects are controlled, 
in our model, by the parameter /oxp (i.e., the rate of ex- 
pansion of Rin after each SN II explosion, see Eq. 1) which 
increases logarithmically with time. It is evident from Fig. [3] 
that, for this particular model, the effect of the expansion is, 
at the beggining, lowering the [Fe/H] content, whereas at a 
later stage the SNe II pollution becomes increasingly domi- 
nant. The [O/Fe] abundance of the forming stars starts in- 
crease monotonically during the evolution, as the Fe content 
is practically constant and SNe II ejects a large amount of 
O (a mean SN II produces 1.82 M© of oxygen). Once the in- 
homogeneous region starts expanding, the O content of the 
forming stars reaches the asymptotic value of [O/Fe]=0.5, 
typical of pure SNe II pollution. 

We use two different values of /exp at two different 
stages, as this gives a better fit to the data. For example, for 
NGC 6752, the initial value of 1.0003 is increased to 1.001 
after ~ 10"^ SNe II have exploded. The effect on the model 
of this variation can be seen in Fig. |3] as an small " hook" 
at [O/Fe]~0.25. Even if it were possible to better tune the 
model (by using another expansion law or simply by finding 
the best parameter combination to maintain [Fe/H] perfectly 
constant during the evolution), we are confident that such 
a fine tuning is not necessary and would not significantly 
change our conclusions. Indeed, given all the other model 
uncertainties we think it is unnecessary to attempt to fine 
tune the expansion caused by the SN II explosions. However, 
to check the robustness of our conclusions to the assumed 
expansion value, we show in Fig. |3] the predicted trend of 
[Fe/H] vs. [O/Fe] for three different values of /exp (roughly 
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Fi gure 7. Th e predicted evolution o f the [Na/Fe] versus [ O/Fe] abundance for NGC 6752, plotted against the observational datasets 
of lYong et alj (2005) (left panel) and ICarretta et al.l |20o3) (right panel, arrows refer to upper limit values; see the original paper for 
details). In our framework the AGB+SN la pollution is due to pre-existing halo stars while only SNe II self-pollute the GC. The solid 
line in the right panel refers to the reference model using the yields of .Karakas fc Latt anzio ( 20071 ) while the dashed lines refer to a model 
with Na production lowered by a factor of four compared with the reference model (see Table [Sj . 
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Fig ure 8. (Left panel) th e predicted evolution of the [Na/Fe] versus [O/Fe] abundance for NGC 2808 against the obser vational data 
from Carretta et al.l 
||2004 ) (squares) and 



2006h. (right panel) the same trend for the case of M 13 plotted against the observational data from lSneden et al.l 



Cohen fc Melendej (diamonds) 



spanning a factor of 10 in dynamic range). It can be seen 
that changing /exp by a factor of three affects the [Fe/H] 
ratio of the forming stars by only 0.1-0.2 dex. 

In Fig. [5] we compare the results from our model with 
NGC 2808 and M 13. These GCs were chosen on the basis 
of their peculiar O depletion with values of [O/Fe] as low as 
— 1.0. These low values are very difficult to obtain with only 
primordial self-enrichment, which may explain the oxygen 
depletion down to a minimum value [0/Fe]~ —0.5 but no 
lower. Extra mixing may also contribute to the depletion 
of oxyg en but only for stars near th e tip of the first giant 
branch. iD'Antona &: Ventural l|2007l ) invoke a combination 
of the two to explain the case of M 13, while in the model 



of iBekki et al] (|2007h such extreme values are achieved by 
assuming that no third dredge-up occurs in intermediate- 
mass AGB stars. 

We are able to reproduce the low values of [O/Fe] with- 
out any particular assumption (see Fig. [SJ. Indeed, such 
low [O/Fe] values arise naturally if the inhomogeneous re- 
gion - polluted by the initial SN la - is smaller in the 
case of NGC 2808 and M13 than NGC 6752 (see Table [T|. 
In the case of M 13, the extreme low values of [O/Fe] of 
— 1.0 requires a slightly smaller Rin as well as less oxy- 
gen from the SNe II that pre-poUuted the halo. To fit the 
properties of M 13, the metallicity of the original ISM gas 
at the time of GC formation needs to be very low, with 
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a typical [ Fe/H] abundance an order of magnitude lower 
([Fe/H]isM = -3.50) than for NGC 2808 or NGC 6752. This 
is despite the final [Fe/H] of M 13 is similar to NGC 6752. 
The lower initial of [Fe/H]isM for M 13 than for NGC 6752 
could indicate that this GC formed at an earlier epoch and 
that it is older than the other two GCs. Photometric ages of 
M M13 compared to jSTGC 6752 suggest that M 13 is older 
l|Rakos fc Schombertllioosi ). but these age determinations 
are highly uncertain and should be treated as tentative. 

In Fig. |6] we show the evolution of the logarithm of the 
metallicity {[og{Z / Zq) versus the [0/Fe] abundance inside 
the inner region for the above models (solid lines) . As in the 
case of the [Fe/H] abundance, the total metallicity remains 
approximately constant during the evolution. The dashed 
lines represent the same evolution without taking into ac- 
count the contribution of metals by intermediate-mass AGB 
stars to the total metallicity fraction. The dotted lines rep- 
resent the case in which the contribution of SN la and AGB 
to the total metallicity are ignored (that is, only SNe II con- 
tribute to the global metallicity). Most of the initial metals 
(e.g., C, N) inside the inner region are due to localized AGB 
pollution, while the SN la contribution is minimal. Indeed, 
even if a SNe la ejects large amounts of Fe (0.74 M©) they 
eject only 1.4 Mq of total metals. Most of the "metal" are in- 
stead produced by AGB stars (an intermediate-mass AGB 
star produces a mean value ~0.08 Mq of metals, most of 
which is N; see Table At the end of the evolution, most 
of the metals are produced by SNe II, where all lines con- 
verge to the same value in Fig. [B] This is also demonstrated 
by the asymptotic value of [O/Fe]=0.5 typical of pure SNe 
H pollution. Since the metallicity of the stars is roughly con- 
stant during the evolution this is also the metallicity of the 
forming stars. 

We do not have a parameter that sets the fraction of 
SNe II that pollute the forming stars. Instead we calculate 
that assu ming the initial mass o f the GC is the same as it 
is today (jPrvor fc Mevlanlll993l ). and assuming a constant 
density and metallicity for the gas in the inner region, the 
fraction of SNe II ejecta that directly enriches the inner 30 
(50) pc (i.e., where the GC star should be forming) turns out 
to be very low, of the order of 0.5-1% (3-6%), depending on 
the details of the model. Indeed, most of the metals produced 
by the SNe II are expelled beyond this limit of ^ 50 pc. 



Table 4. Isotopic ratios for the mean yields of SNe II, SNe 
la, and AGB models. Notations are as in Tables [2] and |3] 





25Mg/24Mg 


26Mg/24Mg 


12C/13C 


SNe II (W&W) 
SNe II (KOB) 
SNe II (C&L) 


7.2 xlO-3 
9.9 xlO-3 
2.8 xlO"^ 


6.3 xlO"^ 
8.9 xlO"^ 
5.0 xlO-^ 


10* 
10^ 
5x10* 


Model (SNe II) 


10~2 




lO* 


AGB (KA07) 


6.5 


20.7 


6.0 



inhomogeneous region expands, the chemical properties of 
the forming stars evolve toward "normal" [Na/Fe]-poor and 
[0/Fe]-rich. In the figure we show two different predictions 
using 1) our reference mod el yields, and 2) the same yi elds 
but with the Na yields from lKarakas fc Lattanziol l|2007h re- 
duced by a factor of four. While the reference model shows 
a clear 0-Na anti-correlation, qualitatively in agreement 
with the observations, it overestimates the initial value of 
[Na/Fe] by ^0.4-0.6 dex. In comparison, the model with 
the reduced Na yields reproduces better the observed anti- 
correlation. Note that the reduced value ((Na)=4.0 x 10"'' 
M q) is in better agreem ent with the Na production reported 
bvlFenner et all (|2004|) ((Na ) = 5.9 x 10"* Mq Mq), and by 
iKarakas fc Lattanziol l|2007l ) for more metal-rich AGB stars 
((Na)=7.7xl0-'* M0 for Z = 0.004). The reduced Na yields 
are also well within model uncertainties (see Tabled. It is 
also possible to obtain a very good match to the absolute val- 
ues of [Na/Fe] by simply reducing (by same factor of 4) the 
amount of AGB pollution compared to the reference model 
(see Table [1} , because the number of AGB stars has little 
effect on O and Fe (see Table However, as discussed in 
Sec. 131 we use AGB nitrogen production as a reference (see 
§ I5.5|l . Thus, in the following we will always assume a mean 
value of {Na)=4.0 x 10~* M© from intermediate-mass AGB 
stars. 

The [O/Fe] versus [Na/Fe] abundance evolution for the 
models tailored for NGC 2808 and M 13 are shown in Fig. |8l 
along with their respective observational datasets. The O- 
Na anti-correlation is well reproduced for these two GCs, 
including the extremely low values of [O/Fe]. 



5.2 The O-Na anti-correlation 

Intermediate-mass AGB st ars do not produce Fe and O 
l|Karakas fc Lattanziol |2007| ) , hence they have a negligible 
contribution to the evolution of these elements. AGB stars 
do, however, produce a considerable amounts of Na. 

In Fig.[7]we show the evolution of [O/Fe] versus [Na/Fe] 
abundances predicte d by our mo del, together with the ob- 
servat ional data from lYong et al.. (2005) and Carretta et al.l 
l|2007l . ; upper limits are shown as arrows) for the case of 
NGC 6752. As already mentioned, the first stars to form 
in our model are [Na/Fe]-rich (due to the Na-rich AGB 
pollution) and [0/Fe]-depleted (due to the localized iron- 
rich SN la ejecta). The effect of the SN la is to mitigate 
the increase of the [Na/Fe] abundance that would other- 
wise reach much higher values. For example, in Fenner et 
al. (2004) the [Na/Fe] abundance increased to [Na/Fe]~1.7. 
Once the first SNe II starts to pollute the ISM and the 



5.3 The Mg-Al anti-correlation 

In the left panel of Fig. [9] we compare the predictions of 
our reference model (solid line) for NG C 6752 with the ob - 
served Al-Mg anti-correlation found bv lYong et al.l (|2005h . 
This dataset shows anti-correlations in stars that are both 
brighter and fainter than the first giant branch bump, 
strongly suggesting that the peculiarities should be present 
in the gas from which the stars formed, and not caused by in- 
ternal stellar evolution. It is apparent that there is a spread 
of ~ 1.3 dex in Al relative to Fe but only a modest spread 
of ~ 0.3 dex in Mg relative to Fe (the corresponding [O/Fe] 
spread discussed in Section [5.21 is ~ 1.0 dex). 

In our framework, the first stars to form have en- 
hanced Al ([Al/Fe]~ 1.4) and depleted Mg ([Mg/Fe]~ 0.2), 
while subsequent generations evolve toward [Al/Fe]~ 0.0 
and [Mg/Fe]~0.5. It is apparent that our reference model 
(Fig. [9)l can only account for the Mg variation and fails to 
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Figure 10. Upper panels: the evolution of the [Al/Fe] versus [Mg/Fe] abundance for NGC 2808 (left panel) and M 13 (right panel). As 
in Fig. |9] the solid lines refer to the model using yields from Karakas & Lattanzioj l|2007t ). whereas the dashed lines refer to the model 
with increased Al from AGB stars. The evolution is plotted against the ob servat ional datasets ofJCarrctta et al. (2006) for the case 
of NGC 2808 (left panel) and ISneden et al.l l|2004l . asterisks^. Johns on et al.l lj2005l . squares) and lCohen &: Melended (|2005l . crosses) for 
the case of M 13. Regarding M 13, we point out that there is an offset of 0.2 dex for the Mg abundance between our model and the 
observational dataset. Bottom panels: the predicted [(Al+Mg)/Fe] versus [O/Fe] abundances; line styles are as above. 



reproduce the ~ 1.5 dex spread in Al. The sm all amount of 
Al produced in the models of Karakas & Latt anzid (|2007l ) is 
not enough to countervail the large amount of Fe deposited 
by the single SN la: the net results is that the [Al/Fe] ratio 
actually decreases in ou r model. Note that even in the case 
of pure AGB pollution (|Fenner et al.|[20o3 ). the maximum 
[Al/Fe] value is underestimated by a factor of ~ 0.7 dex. If 
we assume that low-metallicity AGB stars produce a mean 
value of Al~ 3.5 x 10~^ Mq (instead of the literature values 
of ~ 10"* M0) shown in Table[3l the Al-Mg anti-correlation 
is remarkably well reproduced. 

In the upper panels of Fig. 1101 we plot the [Mg/Fe] ver- 
sus [Al/Fe] evolution of our reference models for M 13 (right 
panel) and NGC 2808 (left panel) against the correspond ing 
obse rvational datasets (Sncdc n ct a l. 2004; Johnson e t al.l 
[2OO5I : ICohen fc Melendej|2005l : ICarretta et al.ll2006l ) . Again 
the reference models are not able to reproduce the Al en- 
hancement, while models with a factor of ~ 50 more Al 
production from intermediate-mass AGB stars are able to 
reproduce the anti-correlation perfectly. The dotted line in 
the right panel represents the same dashed line but offset 
by 0.2 dex in [Mg/Fe]. At low [Al/Fe] values there is a dif- 
ference of ~0.2 dex between the [Mg/Fe] values of the two 
GCs. 



At this point, it is worth asking if our proposed Al yield 
(~ 3.5 X 10~^ Mq) from intermediate-mass AGB stars is 
consistent with the uncertainties present in the theoretical 
models, as well as with other available observations in other 
systems. To test this latter point, we ran a test chemical evo- 
lution model of the Milky Way using GEtool (|Fenner et al.l 
.2002.) with the same enhancement in Al AGB yields that is 
required for the GC model. The result is only a small dif- 
ference (within ~ 0.2 — 0.3 dex and only for [Fe/H]^-1.6) 
compared with the canonical MW model. 

The yields of '^''AI are affected by variations in the 
nuclear reaction rates, as well as the treatment of mass 
loss and convection. For example both, the ^®Mg(p,7)^^Al 
and ■^^Al(p,7)^^Si reactions play an important role in de- 
termining th e final '^'^Al abundance in a given AGB model. 
Ilzzard et al] l|2007t ) concluded that the yields of ^''Al and 
^^Al can vary by up to ~ 2 orders of magnitude, de- 
pending on the mass, metallicity and choice of Mg-Al 
reaction rates. The yields of Al decrease with increas- 
ing mass loss, because more mass is lost before signifi- 
cant Al can be synthesized from Mg via the HBB. The 
yields of nitrogen and Al (but also C, Na and Mg) de- 
pend als o on the third dred ge-up (| Ventura &: D'Antonal 
l2005al lbl: iKarakas et all l2006bl ) . While there is uncertainty 
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Figure 9. Upper panel: the predicted evolution of tlie [Al/Fe] 
versus [Mg/Fe] abund ance fo r NGC 6752. Also shown is the obser- 
vational dataset of Yongjjt aLj t2003). The solid line refers to the 
model with canonical Al yields from iKarakas &: Lattanziol (|2007l ) 
while the dashed line refers to our assumed mean Al production 
{Al=3.5 X 10^ M0) from AGB stars that is required to reproduce 
the observations. The dotted line refers to the previous model but 
with a 50% enhancement of the -^^Mg and ^^Mg yields from AGB 
stars (see Section 15.41 for more details of this model). Bottom 
panel: the predicted evolution of the sum [{Al+Mg)/Fc] versus 
[O/Fe] for the three models described above; line styles are the 



over the amount of dredge-up in intermediate-mass AGB 
stars, observations s uggest that it should be fairly efficient 
l|Wood et all 1 19831 : iGarcfa-Hernandez et al.1 I2OO6I ). More- 
over, hot-bottom burning is the major production site for 
Al (as well as for Na and Mg). In Fig. |5]we show that in- 
creasing the AGB yields of Mg by 50% (particularly ^^Mg 
and 2'^Mg, see next subsection) does not have a significant 
efi^ect on the ability of the model to reproduce the Mg-Al 
anti-correlation observed in NGC 6752. It is not clear, how- 
ever, if the model uncertainties would allow for an increase 
in the Al abundance by a factor of 50 and further study of 
this problem is required. 

Using data for 18 stars, ICarretta et al.l (|2007l 'l noted 
that the sum of Mg-l-Al is approximately constant in 
NGC 2808 for stars in all evolutionary phases (that is, con- 
stant over the whole magnitude range). The conclusion to 
be drawn from this observation is that there has been a 
reshuffiing of Mg into Al in these stars. The bottom panels 
of Fig. |9] and Fig. [10] show the sum of Mg+ Al of our mod- 



els plotted against the [O/Fe] abundance. If we assume the 
above mentioned Al enhancement in AGB yields, the sum 
Al-|-Mg appears to be surprisingly constant (in the reference 
model this sum is only constant within a factor of two) . In- 
creasing the yields of Mg by 50% does not change this result, 
and the sum of Mg+Al is constant to within 0.1 dex (dotted 
line) . 

5.4 The Mg isotope ratios 

A further important test for our model is determining how 
well it can reproduce the Mg isotopic ratios o bserved in 
GC stars (e.g. for NGC 6752 lYong eFal] l2003l ). The Mg 
production from AGB stars is mainly in the form of the 
neutron-rich magnesium isotopes, ^^M g and ^^Mg, which are 
prod uced in the He-burning shell (e.g. iKarakas fc Lattanziol 
l2003h . In Table H we summarize the mean ^^Mg/^^Mg and 
'^^Mg/^'^Mg ratios for the in termediate-mass AGB yields of 
IKarakas fc Lattanziol l|2007l ). as well as the same ratios for 
different SNe II models computed by different authors. It 
is apparent that the ratios ^'^Mg/^^Mg and ^^Mg/^^Mg are 
much higher in the case of AGB pollution compared to SNe 
II, where SNe II preferentially produce ^"^Mg at low metal- 
licities. In the following we will adopt a constant value of 
0.01 for both ^^^Mg/^^Mg and ^'^Mg/^^Mg to reflect SNe II 
(see Table|4| adopting values from different authors does not 
significantly change the results). 

In Fig. [TT] we show the observe d Mg isotopic ratios 
for NGC 6752 from lYong et alj l|2003h . plotted against the 
[O/Fe], [Na/Fe], [Mg/Fe], and [Al/Fe]. In this figure we in- 
clude the predictions from our reference model (solid lines) . 
The agreement is satisfactory, except in the case of Al (solid 
line) which, as already discussed, can be accommodated by 
assuming a larger Al production from AGB stars (dashed 
line). Both the ^^Mg/^'^Mg and ^^Mg/^'^Mg ratios decrease 
as the GC stars evolve from the AGB-SN la peculiar pol- 
lution (2^Mg/^'*Mg~0.5 and ^"Mg/^^Mg^ 0.15), to pre- 
dominantly self pollution by SNe II (^^Mg/^^Mg-O.Ol and 
26Mg/24Mg~ 0.01). 

All GC chemical evolution models so far have failed to 
reproduce the Mg i sotopic ratios observed in NG C 6752 (e.g. 
iFenner et a"i]|2004l : [Ventura fc D'Antona|[2005d ). One of the 
strongest conclusions from lFenner et al. 1 (|2004l i was that it 
is not possible to reproduce the Mg-Al anti-correlation and 
the Mg isotopic ratios with a classical AGB self pollution 
scenario. This is because AGB stars produce Mg instead of 
destroying it, and thus Mg and Al are correlated. In our 
model the [Mg/Fe] depletion is mainly caused by the effect 
of the SN la, hence the production of Mg from AGB stars 
is not significant. The Mg production instead helps to ex- 
plain why strong [O/Fe] depletions (of up to —0.5) are not 
accompanied by a similarly low [Mg/Fe] abundances. Even 
an increase of 50% in the ^^Mg and ^'^Mg yields from AGB 
models does not effect our conclusions, and instead melio- 
rates the Mg isotopic ratios when plotted against [Mg/Fe] 
(see Fig. [TT|) . 

Note that varying the Mg ratios in the SNe II yields 
from 2^Mg/^*Mg and ^^Mg/^-'Mg ~ 0.01 to 0.1 would help 
to fit the observational data. In Fi g. 1121 we sh o w the Mg 
isotopes ratios versus [Na/Fe] from lYong et al.1 (|2006al ) in 
the case of M 13 together with the predictions from our 
model and, again, the observations are matched fairly well. 



The Chemical Evolution of Globular Clusters 13 



1.0 

0-8 
^ 0.6 
0.4 

CO 

" 0.2 
0.0 



0.4 
^ 0.3 
^ 0.2 

60 

^ 0.1 





1 1 


□ 1 


□ 1 ■ 

II 










. 











-0.5 0.0 0.5 -0.20.0 0.2 0.4 0.6 0.2 0.3 0.4 0.5 0.0 0.5 1.0 1.5 

[0/Fe] [Na/Fe] [Mg/Fe] [Al/Fe] 



0.0 



NGC 6752 




□ 




□ 
□ 



-0.5 0.0 0.5 -0.20.0 0.2 0.4 0.6 0.2 0.3 0.4 0.5 
[0/Fe] [Na/Fe] [Mg/Fe] 



0.0 0.5 1.0 1.5 

[Al/Fe] 



Figure 11. Predicted evolution of the Mg isotopic ratios versus [0/Fe], [Na/Fe], [Mg/Fe], and [Al/Fe] for our model (dashed lines). The 
dashed lines refers to our assumed "Model" (see Table |3j, while the dotted lines represent the evolution in the case of the ^^Mg and 
26 Mg AGB yields increased by 50%. (see Table gjl. 
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Figure 12. Predicted evolution of Mg isotopic ratios versus 
[Na/Fe] for our referen ce model (solid line s) compared with the 
observational data from lYong et al] l|2006al ) . The dashed lines cor- 
respond to the model with 25Mg/24Mg=0.1 and 26Mg/24Mg=0.1 
from SNe II (instead of 0.01; see Tablegll. 



Assuming 25Mg/2*Mg=0.1 and ^^Mg/^^Mg^O.l for the SNe 
II yields does an even better job at reproducing the obser- 
vations. 



5.5 The C-N anti-correlation 

In Fig. [13] we illustrate the evolution of [C/Fe] versus 
[N/Fe] of our reference mo del for NGC 6 7 52, t ogether 
with observational data by ICarretta et al.l |200^ for 9 
dwarfs (filled diamonds), 9 sub-giants stars (open dia- 
monds), a long with 2 giants st ars (open triangles) ob- 
served by ISmith fc NorrisI (|l993l ). Comparison with stars 
near the main sequence is preferred, because variations due 
to the first dredge-up and deep mixing can affect the abun- 
dances of evolved stars. However, it is very difficult to ob- 
tain high-quality spectra of main-sequence stars owing to 
their low brightness. In the following we will not discuss 



Denissenkov fc Toutll2000l : IWeiss et al.ll2000l : ICratton et al l 
and references therein). 
While the more evolved stars (open symbols) show a 
clear C-N anti-correlation, similar to that observed in other 
GCs (see also Fig. [14] for M 13), the inclusion of the dwarf 
stars m akes it difficult to distinguish a clear evolutionary 
pattern. ICarretta et al.l |2005) pointed out that the spectra 
for these turn-off stars did not have enough quality to allow 
for an accurate determination of the carbon abundances; 
the same also applies to the ^^C/^'^C ratios (see 15. 7p . Typ- 
ical errors for the [C/Fe] abundances are of the order of 
~ 0.2 dex with the exception of the dwarf stars that show 
larger uncertainties. 

A quick examination of Table [3] shows that nitro- 
gen is the element produced in the largest proportions by 
intermediate-mass AGB stars. Thus, it is not surprising that 
the first stars forming in the inner region initially polluted by 
AGB stars show very high values of [N/Fe] ( ^ 1.5) despite 
the action of the SN la. For comparison, in iFenner et al.l 
(2004), stars polluted with AGB ejecta show [N/Fe] values 
as high as ~2.0). Carbon is also produced in intermediate- 
mass AGB stars, although not to the same level as N, hence 
the [C/Fe] ratios remains low. A glance at Fig.[l3]shows that 
while a large variation of [N/Fe] is achieved in our model in 
accordance with observations, the ratio of [C/Fe] remains 
almost constant. While [C/Fe]~ 0.0 at high values of [N/Fe] 
is consistent with the observed [C/Fe] abundances inferred 
from dwarf stars, assuming a factor of 4 less carbon is pro- 
duced by AGB stars does a much better job of reproducing 
the C-N anti-correlation observed in sub-giant stars, (see 
Fig. [13] and Table [3}. 

The C-N anti-correlation is also seen for stars on the 
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Figure 13. Upper Panel; evolution of [C/Fe] versus [N/Fe] for 
the reference model (solid line) for NG C 6752, plotted a gainst 
observational data of two AGB stars by ISmith &: Norrig lll993l ) 
(triangle), and 9 sub-giant s stars (open d i amon ds) and 7 dwarf 
stars (filled diamonds) by ICarretta et al] ||2005|) (squares). The 
dashed line represents the model with AGB C yields decreased by 
a factor of 4. Lower panel; evolution of the sum of [(C+N+0)/Fe] 
(solid lines) for the reference model (thick lines) and well as for 
the model with r educed C (thin l ines) p lotted against the observa- 
tional dataset of lCarretta et al] | |2005| ) ; the symbols are the same 
as in the upper panel. The evolution of each single element is also 
shown; C (dotted lines), N (dot dashed line), and O (dashed lino). 

mam sequence turn off of M 13 (jBrilev et all l2004bl ) (see 
Fig. I14|l , whereas we are not aware of a similar dataset for 
NGC 2808 (although the C-N evolution of this GC is simi- 
lar to M 13 for more evolved stars). We refer to lBrilev et alJ 
l|2004bl ) for a discussion of the uncertainties affecting the 
abundances; these can be large (up to 0.5 dex, see their 
Fig. 5) owing to the faintness of the observed stars. Again 
the two models show quite similar results, with [C/Fe] re- 
maining approximately constant in the reference model, and 
the model that assumes C is under produced by a factor of 4 
doing the best job of reproducing the C-N anti-correlation. 



To test the possible contribution of low-metallicity 
fast rotators or mass lo ss in metal-poor massive stars 
(jMevnet fc Maedeij [2002bl lah we qualitatively increase the 
nitrogen mean yield for SNe II, from our reference value of 
[N/Fe]=-0.90 to values of [N/Fe]=-0.30 and 0.0. This in- 
crease only slightly changes the results, by increasing the 
minimum value of [N/Fe]~ —0.4 in Fig. [14] by a factor of 
0.20-0.40 dex. 

5.6 The sum of C+N+O 

One of the most important observational constraints that 
any GC chemical evolution model nee ds to fulfill is the 
constancy of the s um C-hN-hO (e.g llvans et alJ Il999l : 
ICarretta et alj |2005| ) . Self pollution scenarios using AGB 
models calc ulated using the mixing-length theory (MLT) of 
convection (|Fenner et al.t2 004: Kar akas et al. 2006a ) predict 
that stars forming from different generations of AGB ejecta 
show large increases in the CNO sum. This is in contrast to 
AGB models computed using the "full spectrum of turbu- 
lence" convective model (| Ventura fc D'Anto na 2005b) that 
manage to keep the sum of CNO constant within a fac tor of 
two. This is because the IVentura fc D'Antonal (|2005bh show 
little third dredge-up (which increases C and hence N) and 
more efficient HBB (which converts C and O to N). Hence, 
MLT models create C and N rich AGB e jecta at almost con- 
stant O (see Fig. 4 in lFenner et aLll2004l ). FST models result 
in moderately O-d epleted AGB ejecta that is on ly moder- 
ately enriched in N IVentura fc D'Antonal l|2005bl ') 

In the bottom panel of Fig. [T3] we show the C-|~N-|-0 
sum (solid line) of our model for NGC 6752, together with 
[C/Fe] (dotted lines), [N/Fe] (dot-dashed line) and [O/Fe] 
(dashed line) plotted against [N/Fe]. In t he same figure 
we also include the observational data from ICarretta et al.l 
(2005) for 5 sub-giants and 7 dwarf stars from the same GC. 
It is quite remarkable that our model predicts a practically 
constant CNO sum and is a satisfactory fit to the obser- 
vations. It is apparent that while the [C/Fe] abundance re- 
mains approximately constant, a large variation of [N/Fe] 
comes with a similarly large variation in [O/Fe]. The N de- 
posited by AGB stars is counteracted by the O production 
in SNe II when the inner region expands and mixes (lower- 
ing the [N/Fe] value]). In principle, in our model, there is no 
reason for this to happen because the mechanisms control- 
ling the Fe abundance and the CNO sum are different (i.e., 
SNe II, SN la, and AGB). For example, more AGB pollu- 
tion (or no AGB pollution at all) would increase (decrease) 
the initial CNO sum and the condition C-|-N-|-0=constant 
would not be fulfilled. However, in the three particular GC 
considered here, once the free parameters are fixed to re- 
produce the observed anti-correlations, a constant C+N+O 
naturally arises. 

The thin solid lines and dotted lines in the lower pan- 
els of Fig. [TS] and Fig. [T?] refer to the case in which the 
C production from AGB stars is lowered by a factor of 4 
to better fit the C-N anti-correlation. In this case the sum 
of C+N+O also remains quite constant, in agreement with 
the observations. The fact that a variation of 0.4 dex of the 
[C /Fe] abundance only slightly changes the CNO sum is un- 
derstandable if we consider that O and N are a factor of ~ 
8 more abundant than C at low values of [N/Fe] and high 
values of [N/Fe], respectively. 
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Figure 14. M 13 Upper Panel: evolution of [C/Fe] versus [N/Fe] 
for tlie reference model (solid line) p lotted against the obser- 
vational data of iBrilev et al.l l l2004af ). The dashed line corre- 
sponds to the model where we decrease the C yields from AGB 
stars by a factor of 4. Lower panel; the evolution of the sum of 
[(C-|-N-|-0)/Fe] (solid lines) for the reference model (thick lines) 
as well as for the model with reduced C (thin lines). The evo- 
lution of each single element is also shown; C (dotted lines), N 
(dot dashed line) and O (dashed line). The corresponding two 
diagrams for the case of NGC 2808 are almost identical to the 
results for M 13 and not shown. 

5.7 The carbon isotope ratio 

In Fig.[T5]we show the evolution of the ^^C/^^^C ratio versus 
[Na/ Fe] from our model p lotted against observational data 
from lCarretta et al.l l|2005l 'l for sub-giant stars in NGC 6752. 
Owing to the poor quality of the spectra, these authors were 
not able to derive the carbon isotopic ratio for dwarf stars. 
Our model over-predicts the ^^C/^'^C ratio at any value of 
[Na/F e] an d by up to a f actor of 1.5 dex. IChiappini et al.l 
i|2008h and lHughes et all (|2008l 'l note that a model without 
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Figure 15. Carbon isotopic ratio i^C/^'^C versus [Na/Fe] for 
the NGC 6752 reference mod el, plotted against the observational 
dataset of lYong et al. I 1I2OO3I) . The dashed and dot-dashed lines 
represent a model in which the possible contribution of low metal- 
lic ity, fast rotators are qu alitatively taken into account according 
to lChiappini et al.l l|2008l ). 



fast rotators gives very high values of the C/ C ratio of 
up to 10* because SNe II mainly produce '^'^C (solid line 
in Fig. I15|l . Such high values are not reached in Fig. [15] 
because of the partial contribution from intermediate-mass 
AGB stars, that also produce low_^G/^^C ratio yields (see 
Table m. The fChiappini et al.! (j2(joi ) model with rapidly ro- 
tating massive stars brings the ratio down to ^■^C/^''C~80 
(dashed line) at [Fe/H]=-3.5 and down to "C/^^C~30 
(dot-dashed line) at [Fe/H]=— 5.0. If low-metallicity, rapidly 
rotating massive stars contribute some ^"^G, the predicted 
carbon isotope ratio will be lower but not as low as required 
to match the observational data (see dashed and dot-dashed 
line in Fig. [TSjl . 

Evidently it is not possible to reproduce the observed 
12q^13(-i j-g^^jQg without invoking internal stellar evolution 
(e.g., the first dredge-up and/or extra-mixing processes). In 
the more luminous sub-giant stars the observed isotopic ra- 
tio ^■^C/^'^G is very low, typically in the range of ~ 3 to 
11 (noting that the equilibrium value of the CNO cycle 
is ~ 3). The first dredge-up can reduce the ratio from 90 
(solar) to ^ 20 in standard s t ellar evolution m odels (e.g., 
iBoothrovd fc Sackmannlll999t ). IShetrond l|2003h found that 
in NGG 6528 and M 4 this ratio declines steeply with in- 
creasing luminosity along the RGB. The discontinuity in the 
i^G/^G ratio occurs at the bump luminosit y and has also 
been detected in metal-poor field halo stars (jGratton et al.l 
|2004 ). 

iGarretta et al. I l|2005l ) inferred a lower limit of 10 for 
the ^^G/"G ratio for the three dwarf stars in NGG 104. 
Even if these values are more in accordance with our model, 
further determinations of the ^^G/^^G ratio in turn-off stars 
is needed before making a more detailed comparison. 
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Figure 16. Evolution of the helium content, Y, versus [O/Fc]. 
The red solid line represents the model for NGC 6752; the blue 
dotted line for M 13; the dot dashed orange line for NGC 2808. 



5.8 The Helium Content 

The unusual HB morphology observed in NGC 2808, which 
exhibits an extended blue tail and a Ra p separating the red 
and the blue clumps (|Bedin et al.ll200ol ). can be reproduced 
if the blue stars have a higher helium content {Y ~ 0.32) 
compared to t hose in the red clump, wi th primordial helium 
{Y ~ 0.24) (|D'Antona fc Caloil [20o3 '). N GC 2808 is also 
know n to have a peculiar main sequence (iD'Antona et al.l 
I2OO5I ). in which the bluer stars are inferred to have an higher 
helium content from fitting theoretical isochrones to the ob- 
served d ata. To reprodu c e the colour-magnitude diagram 
(CMD), ID'Antona et al.l (|2005l ') suggested that ~ 20% of 
the stars have a helium content as high as F ~ 0.40, ~ 30% 
have a spread of Y between 0.24 and 0.29, while the remain- 
ing 50% of the stellar population have primordial Y . Further 
observations discov ered a triple main sequence in NGC 2808 
iPiotto et al.ll2005l ). that has given support to the extreme 
value of y ~ 0.40 required to fit the CMD. 

Even if such large Y values are required to reproduce 
the main features of the CMD, the origin of the helium is 
not known. Helium values as high as y = 0.40 are very dif- 
ficult to reproduce with canonical chemical evolution mod- 
els without violating other observationa l constraints (e.g. 
iKarakas" et al-lbOOeaT ). ID'Antona fc Caloil { 20041 ) required a 
factor of ~ 10 more intermediate-mass AGE stars than 
what it is inferred from a Salpeter-like IMF to return the 
amount of Y needed to form the number of blue HB stars in 
NGC 2808. Karakas et al. ( 2006a) highhghts the difiiculty in 
producing the large postulated helium enrichment using the 
AGB self-pollution scenario. In particular, their largest pre- 
dicted value of y ~ 0.29 is accompanied by a large increase 
in sum C-l-N-l-0. 

In Section FS.GI we show that one of the main successes of 
this model is the ability to reproduce the C-l-N-l-O sum for 
all three GCs NGC 6752, M 13, and NGC 2808. In Fig. [TS] 
we plot the helium content Y of the reference models for 
each cluster against their [O/Fe] abundances. Note that we 



assume a primordial value of y ~ 0.248 (|Coc et al.ll20o3 ). It 
is quite evident that in the cases of NGC 6752 and M 13 the 
increase is quite modest and within Ay ~ 0.015, whereas 
for the specific case of NGC 2808 we obtained Y ~ 0.29. 
This larger rise in the Y values is due to the larger density of 
polluting intermediate-mass AGB stars in NGC 2808, which 
is ~2-3 times greater than for the other GCs. 

Even if NGC 2808 is the only GC out of the three that 
shows a appreciable helium variation, a maximum value of 
y ~ 0.29 is smaller than the postulated very high value of 
y ~ 0.40 required to e xplain the most en riched stars in the 
triple main sequence l|Piotto et al.l |2007|). It is very c lose, 
however, to the value of y ~ 0.32 (" Bedin et al] |200G| ) re- 
quired to explain the unusual HB morphology for NGC 2808. 



6 COMPARISON WITH PREVIOUS MODELS 
AND GENERAL COMMENTS 

To explain the peculiar chemical properties observed in GCs, 
we propose a model based on a new idea of pre-peculiar in- 
homogeneous pollution produced by the explosion of a single 
SNe la. In principle, due to the low star formation rate dur- 
ing the formation of the proto-halo, the possibility that 2 or 
3 (or more) SNe la explode essentially simultaneously (<20- 
30 Myr apart in time) in a region smaller than 30-50 pc is 
much less probable than having a single event. Having said 
that, our model can accomodate such a scenario without a 
dramatic impact on the predictions based upon the base sce- 
nario of a single SNe la event. For example, if 2 (or 3) SNe 
la explode in the same region, ejecting therefore a factor of 
2 (or 3) more iron, an expansion (i.e., an increase of Rin) 
of a factor of ^ (or v^) would lead to the same [Fe/H] 
and chemical properties seen in the base model. In the case 
of no SNe la, but intermediate-mass stellar AGB pollution, 
one is lef t with the classica l "external" GC self enrichment 
scenario (|Bekki et al.ll2007l ) in which the high values of N, 
Na (and other elements produced in substantial amounts by 
AGB stars) are easily explained, but it is very difficult to 
explain the depletion in oxygen. Since all the well-observed 
GCs seem to have an oxygen-depleted population of stars at 
some level, in our framework we must assume that the inho- 
mogenous contribution of SN la is the seed for the formation 
of GCs. 

Of course, a larger number of SNe la may contribute to 
the enrichment of the proto-halo ISM (outside the proto-GG 
region). We believe, though, that the approximation of an 
ISM being enriched only by SNII products is sound one, as 
the SNe la contribution to the "mean" MW's stars' chem- 
ical properties starts being important only after ~1.5 Gyr 
and at [Fe/H] > -1 (Tinsley 1979; Edvasrdsson et al. 1993; 
McWilliam 1997). The [Fe/H]isM values in our models are, 
in a ll cases, lowe r than this value. 

iBekki et al.1 (|2007h point out that one of the greatest 
advantages of an external AGB-poUution scenario is that 
the fraction of peculiar stars observed in GCs is not limited 
to the amount of a stars forming in the first generation. In 
particular, in this model the fraction of O-depleted to O- 
rich stars depends only in the way that the SNe II pollution 
proceeds, while the number of AGB stars initially polluting 
the inner region is quite low (JVagb ^ 200; see TablelU. The 
total amount of stars needed to create a typical value of 200 
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intermediate-mass AGB stars is ~ 10'' Mq , which is a smaU 
fraction of the total MW Halo or GC mass. Even assuming 
that this number of stars are trapped in the potential well of 
the forming GC and are homogeneously distributed within 
50 (or 100) pc, their density would be 1.9 x 10"^ (2.4 x lO'^*) 
M0pc~^ which is low compared t o the densities observed in 
GCs (e.g.. |Prvor fc Meylan|[l993h . 

In principle, our proposed scenario has no particular 
problems explaining the observed high fract ion of peculiar 
to no rmal stars observed in some GCs (e.g.. [Carretta et al.l 
I2OO6I . for NGC 2808). The self-pollution scenario, on the 
other hand, requires ad hoc assumpti ons such as an a,noma- 
lous IMF peaked at 4--8 Ma (e.g.. ISmith fc NorrisI Il982l : 
iD'Antona fc Caloil [20oi : iPrantzos fc Charbonnell 120061 ): or 



that the GC was initially much more massive and has lost 
~90-99 per cent of the stars that were formed in the first 
generation through tidal stripping. Our model does not re- 
quire these assumptions: the only assumption we require is 
that star formation should proceed more or less smoothly 
while SNe II explode and self enrich the proto-GC; and the 
chemical properties of the forming stars evolve from peculiar 
to normal. The fraction of normal to peculiar stars depends 
on how the SF proceeded and we will try to address this 
problem with detailed hydrodynamical simulations. 

In previo us rno dels (e.g. Fenner et al.l |2004| : 

iD'Antona eFal] l2005l : iBekki et all boOTi ) the CN bi- 



modality is not reproduced due to the lack of C-depleted 
and N-rich stars, since C is produced in AGB stars. The 
same is also true for Mg, with the result that positive C-N 
and Mg-Al correlations a re found in contra diction to the 
observational constraints. iBekki et al.l (|2007l ) tried to solve 
this problem by using yields from AGB models that had no 
third dredge-up (as these authors noted, this assumption is 
purely hypothetical). The smaller production of C is what 
allows this study to reproduce the C-N anti-correlation. 
The O-Na and Al-Mg trends were not recovered because 
the lack of dredge-up resulted in a smaller production of Na 
and Mg (mainly ^®Mg). Note that in all AGB models, O 
destruction via HBB is not efficient enough to account for 
the very low O abundances observed in some GC stars (typ- 
ically [Q/Fe] ~ —0.5 or less ). The no n-standard models of 
Ventura fc D'Antonal (|2005d ) (see also lVentura fc D'Antonal 



2006h that use a different convection model, can reproduce 



most of the abundance anomalies but only when including 
convective overshooting to force some third dredge-up along 
with adopting high mass- loss rates. In our model none of 
th e above assumptions a r e requ ired. We use the AGB yields 
of iKarakas fc Lattanziol l|2007l ) , averaged over a canonical 
Salpeter IMF. 

The success of our model is the combined effect of in- 
homogeneous AGB pollution together with the effect of a 
single SN la and its Fe-rich ejecta. The net effect is to lower 
the ratio [X^/Fe], where the element Xi is negligibly pro- 
duced (or destroyed) relative to SN II in intermediate-mass 
AGB stars (e.g., oxygen). In the case where Xi is produced 
by AGB stars, the ratio is either maintained at an approxi- 
mately constant value or the enhancement is mitigated ow- 
ing to the effect of the SN pollution (e.g., C, Na, Mg and 
possibly Al). The only real problem that our model suffers is 
the production of Al should be enhanc ed by a factor of ^ 10- 
50 compared to the AGB yields from lKarakas fc Lattanziol 



(|2007l ). Our model, however, is the first that has shown any 
promise at explaining the Mg isotopic ratios in GC stars. 

One possibility in a canonical AGB self-pollution model 
is that the second generation of stars did not form from 
AGB-enriched material but t hese stars had their surfaces 
polluted by AGB winds (e.g.. [Pantona et al.lll983l ). In our 
framework the O-depleted stars obtain their abundances 
from SNe la pollution mixed with AGB material. We could 
image a similar surface pollution model as follows: after the 
first "normal" generation of stars have formed, self pollu- 
tion by intermediate- mass AGB stars and SNe la (which has 
similar timescales; see Section[21) start inhomogeneously pol- 
luting the surface of stars. It would be difficult, however, to 
explain all the chemical properties of GC stars with surface 
pollution. For example, it would probably not be possible 
to maintain [Fe/H] constant in the most O-depleted stars. 
While surface pollution of SN la ejecta would decrease the 
[0/Fe] ratio it would also dramatically increase the surface 
Fe content of the star. 



7 CONCLUSIONS 

We have modelled the chemical evolution of GCs under the 
hypothesis that the abundance anomalies were set during 
the formation of the Milky Way Halo (or inside their host 
dwarf galaxies) . Inside a volume of gas enriched to some level 
by a first generation of low-metallicity SNe II, GC forma- 
tion takes place owing to the inhomogeneous (respect to the 
surrounding ISM) effect of a single SN la plus material from 
intermediate-mass AGB stars. After a burst of star forma- 
tion, new SNe II start to explode, self-polluting and expand- 
ing out the inhomogeneous region, which mixes with the 
surrounding (metal poorer) ISM. In our framework, all the 
peculiar stars observed in GCs (e.g., O-depleted, Na-rich) 
were born in the above inhomogeneously enriched volume, 
while "normal" stars (i.e., 0-rich, Na-depleted) are formed 
successively once the SNe II associated with the SF have 
washed out the inhomogeneous region and self-polluted the 
ISM. In conclusion, our main findings can be summarized as 
follows: 

(i) We reproduce the O-Na anti-correlation. In our model 
the first stars to form have low [O/Fe] and high [Na/Fe] 
abundances, with the system then evolving toward "normal" 
abundances once SNe II start polluting the ISM. 

(ii) We reproduce the Mg-Al anti-correlation, but 
only by assuming that intermediate-mass AGB produce 
more Al than predicted by theoretical models (e.g., 
iKarakas fc Lattanziol |2007| ) . This assumption needs to be 
carefully studied in the context of AGB model uncertain- 
ties. 

(iii) We satisfactorily produce the Mg isotopic ratios ob- 
served in NGC 6752. 

(iv) The [Fe/H] abundance remains approximately con- 
stant during the evolution. This is because of mixing caused 
by the subsequent SNe II explosions between the Fe-enriched 
inhomogeneous inner region and the [Fe/H]-poorer ISM. The 
Fe produced by the SNe II themselves does not significantly 
affect the resulting Fe abundance in the forming GC stars. 

(v) We reproduce the C-N anti-correlation within obser- 
vational uncertainties, and the sum C-l-N-l-O remains ~ con- 
stant, in agreement with the observations. The best fit is ob- 
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tained, however, by assuming that the AGB stars produce a 
factor of four less C. This assumption may not be in accord 
with current theoretical models. 

(vi) The model cannot reproduce the ^^C/^^C ratios ob- 
served in evolved giant stars in the GCs, and internal stellar 
evolution processes are required. 

(vii) GCs are a factor of ~ ten more metal rich than co- 
eval MW halo stars (which are the assumed stellar host pop- 
ulation). 

In a forthcoming paper we will focus on other elements 
(eg, Li, F, Si, Ca, Ti, V, Co, Ni, and Cu) to further tost our 
framework. We will also focus upon the extreme metallicity 
GC cases that may prove to be a challenge for our model: the 
very metal poor M15 ([Fe/H]=-2.26) and the metal-rich 47 
Tuc ([Fe/H] =— 0.76). Preliminary results are encouraging. 

One last question wc could ask is why would all GCs 
show such a peculiar enrichment as outlined in this paper? 
Even if we do not have a conclusive answer, the region en- 
riched by the SN la and AGB stars has a [Fc/H] abundance 
and metallicity, Z, that is a factor of ~ ten higher than the 
surrounding ISM. Thus could be due to the corresponding 
enhanced cooling, this region may be more likely than the 
halo field to produce a star cluster with regions of enhanced 
star formation. In our framework, the SN la (and AGB) pol- 
lution is the seed which may have led to the formation of 
the globular cluster itself. 
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